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Abstract. A standard assumption in quantum chaology is the absence of correlation 
between spectra pertaining to different symmetries. Doubts were raised about this 
statement for several reasons, in particular, because in semiclassics spectra of different 
symmetry are expressed in terms of the same set of periodic orbits. We reexamine 
this question and find absence of correlation in the universal regime. In the case of 
continuous symmetry the problem is reduced to parametric correlation, and we expect 
correlations to be present up to a certain time which is essentially classical but larger 
than the ballistic time. 
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1. Introduction 

The topic we are to consider is the cross-correlation between spectra of different 
symmetry for so-called quantum chaotic systems, i.e. quantum systems, whose classical 
analogue is chaotic. The instinctive reaction would be to say that states of different 
symmetry live in different spaces, don't talk to each other and cannot conceivably 
display any correlation |T] . Doubts appear when one remembers that in the semiclassical 
approach the level correlations of different symmetries can be expressed in terms of sums 
over the same set of periodic orbits, and such "brotherhood in the parent orbits" may 
introduce a correlation. Although independence is habitually assumed, the question 
deserves a closer look. In nuclear physics some evidence has been assembled that 
correlations in two-body random ensembles with finite sets of states exist, which is not 
too surprising considering the small number of independent two-body matrix elements 
[2], and certain nuclear scattering data support the existence of correlations to some 
extent |3] but a study of correlations between spectra of opposite parity in the limit of 
many particles revealed no correlations 

These findings make it all the more desirable to analyze the usual assumption of 
independence in the semi-classical limit in some detail, as we shall do here. We get 
basically a negative result, both for the case of discrete and the continuous symmetry, 
i.e., the absence of correlations is confirmed in the universal regime. 

These conclusions are to be taken with a grain of salt. Correlations may still 
exist in certain pathological systems (mostly those with disconnected phase space often 
due to unusual behavior under time- reversal symmetry |5l El [7]. Our reasoning also 
does not necessarily apply to arithmetic billiards which have chaotic dynamics, but an 
exponentially large number of geodesies with identical actions [8]). Nevertheless, it has 
been shown by explicit evaluation of the Selberg trace sum formula for the correlation 
function, that eigenvalues for even and odd parity in the modular domain on the surface 
with constant negative curvature are indeed uncorrelated [9] More importantly, for 
continuous symmetries correlations may exist up to a certain time which is essentially 
classical in nature, if the difference of quantum numbers characterizing the irreducible 
representations is of order one as the quantum numbers themselves get large. Note that 
this effect can be easily overlooked, as times are often given in terms of the Heisenberg 
time which, in the semi-classical limit, goes to infinity. 

The paper is structured as follows: After laying foundations we discuss the discrete 
symmetries in what we consider to be the simplest case. This is the reflection symmetry 
in a two-dimensional system when we can pass to a half-space with Dirichlet or Neumann 
boundary conditions. We obtain the relevant results for these two particular cases 
using periodic orbit expansions and find the expected absence of correlations. Starting 
with the diagonal approximation we discuss later the off-diagonal contributions and 
times beyond the Heisenberg time. After that we study general discrete symmetries. 
Trivial generalization of the described technique is possible only for groups that induce a 
decomposition of space into fundamental domains. This is not always true and we shall 
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present the general argument using the symmetry decomposition of periodic orbits, as 
proposed in [TO] . 

In the next section we pass to continuous symmetry and consider the effect of a 
change of the index of irreducible representation as parametric correlation [11] within the 
symmetry-reduced Hilbert space. Again, we find absence of correlations in the universal 
regime, i.e., on the scale of the Heisenberg time, however correlations on a classical i.e. 
/^-independent time scale are expected. Finally we present some conclusions and an 
outlook about possible observations. 



2. Definitions: spectral cross-correlation function and the form factor 

We start with some definitions that will be essential to the later arguments, and at the 
same time we fix the notations we shall use. 

Let p^^\E) = and p^^\E) = Efc^(^-^f^) be two spectral 

densities describing either two subspectra of different symmetry of the same Hamiltonian 
or the spectra of two different Hamiltonians. We shall assume that both spectra are 
unfolded to the same constant mean level density A (see, e.g., [1]). The precise manner 
of the unfolding procedure is a somewhat thorny issue and will be discussed later. 

The two-level cross- correlation function is defined as, 

" — mEwmm — ' 

where (...) denotes averaging over an interval of E and smoothing over a window of 
the dimensionless energy offset e; since the spectra are unfolded we have (p^^-* (E)) = 
{E))=A-\ 

It is often more convenient to work with the Fourier transform of R'^^''^\E,e) with 
respect to e. The result is the cross form factor which is a double sum, 

ir(i'2)(r) = 5^e<^^''-^''>^-/^ (2) 

k,l 

the dimensionless r = T/Th is time in units of the Heisenberg time Th = 27r^/A. 

A semiclassical representation for the spectral correlator and the form factor follows 
from the Gutzwiller expansion of the fluctuating part of the spectral density in terms 
of the classical periodic orbits 7, 

Posc(^)~5^Ae*''^/''-*'^--/2_ (3) 

7(1) 

Here S^, A^, p^ are the action, stability coefficient and the Maslov index of the orbit 7. 
Strictly speaking this expression should also include the periodic orbit repetitions; these 
however become irrelevant in the semiclassical limit. Replacing in ([2]) both spectral 
densities by the respective expansions ([3]) we get, 

K(i,2) (r) ~ 1 / Yl A,,A7,e<"- )/^'-<^- -^->/A ; (4) 

\Tht<T~,, , <Tht+S I 
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Figure 1. Cardioid billiard. Fundamental domain is bordered by the symmetry line 
(dashed) and the right (bold) half of the billiard wall. 

the brackets (...) denote averaging over E and smoothing over a small window of r. 
The sum is taken over periodic orbits 71 of the system 1 and 72 of the system 2. Their 
periods , must lie in the interval [T, T + 5] whose width 6 is small compared with 
T = tTh. 

3. Discrete symmetry 

To illustrate the central point of our argument we shall start with the simplest possible 
example of a system with a single reflection symmetry. 

3.1. A simple example: reflection symmetry 

Consider the spectral correlation between the even and odd spectra in a chaotic system 
with two degrees of freedom such as the cardioid billiard [I2l [131 Ej, see Fig. [H On the 
line of symmetry the even and odd wave functions either have zero normal derivative 
or themselves turn into zero. The subspectra of definite parity can be found separately 
by solving the Schrodinger equation in the desymmetrized billiard (the fundamental 
domain) obtained by cutting the billiard along the symmetry line and dropping one of 
the halves. Imposing the Neumann (resp. Dirichlet) boundary condition on the cut we 
shall obtain the even (resp. odd) part of the spectrum. 

The semiclassical level density consists of the smooth part (the Weyl term) and the 
fluctuating part which is the Gutzwiller sum over periodic orbits of the desymmetrized 
billiard with the mirror- reflecting wall at the symmetry line. The difference between the 
even and odd cases in the Weyl term is responsible for the fact that the A^-th eigenvalue 
in the symmetric sector Efj lies lower than the A^-th eigenvalue in the antisymmetric 

/—A S 

sector by an amount of the order \/Ej^ , see [15] and reference therein. This issue 
has some consequences but can be taken care by an appropriate unfolding procedure. In 
the Gutzwiller sum the even and odd cases differ by the effective Maslov indices of the 
periodic orbits since each reflection in the Dirichlet and Neumann boundary contributes 
71 or zero, respectively, to the phase of the orbit contribution. 
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We shall denote by the number of strikes of the periodic orbit 7 of the 
desymmetrized billiard against the line of symmetry. It is easy to see that orbits with 
even correspond to the periodic orbits of the full billiard with the same period; 
the latter are non-symmetric unless 7 is self-retracing. On the other hand orbits with 
odd are associated with symmetric orbits of the full billiard whose period is doubled 
compared with 7. It is physically reasonable to assume (although we don't know whether 
the result is proven) that the number of orbits 7 with periods in a certain interval 
[T, T + 6] with T large may not depend on parity of n^, as well as on parity of the 
number of visits to any side of the desymmetrized billiard. 

The semiclassical cross form factor now becomes a double sum over the same set 
of periodic orbits of the desymmetrized system. In the following expression summation 
over 71 and 72 is connected with the Gutzwiller expansion of the even and odd spectral 
densities, respectively, 



\Tht<T^^,T-,^<Tht+5 I 

the factor (— 1)"^2 emerges because each visit to the symmetry line with the Dirichlet 
boundary condition yields the Maslov phase vr. 

3.2. Diagonal approximation 

Dropping in ^ all summands with 71 not coinciding with 72 (up to time reversal since 
we have time reversal invariance) we come to the diagonal approximation. 



It differs from Berry's diagonal approximation for autocorrelation form factor [T6] by 
the presence of the sign factor (—1)''^ leading to the destructive interference among the 
contributions. Assuming that for the long orbits with very large the number of orbits 
with even and odd is the same and not correlated with the stability coefficients we 
obtain i^g„,diag (^) = 0. 

A warning has to be made here. Our basic argument (even and odd spectral 
densities are created by the same set of orbits which contribute with the same action but 
a different Maslov index) implicitly assumes ergodicity of the classical motion and loses 
its validity if ergodicity is violated. Consider, e.g., a particle moving in a symmetric 
double-well potential in two dimensions. For energies below the dividing barrier the 
classical motion is clearly non-ergodic since the available phase space disconnects into 
two symmetric domains; even and odd levels in the quantum problem are then almost 
degenerate, up to the tunnelling splitting. Consequently the cross g < — > u form factor 
will practically coincide with the autocorrelation form factor of the spectrum in a single 
well. 
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Another trivial issue should also be pointed out here: for very short times (ballistic 
times) corresponding to the period of the shortest classical orbits, there are no more 
cancellations between different orbits and correlations certainly arise. 

3.3. Off-diagonal contributions 

The off-diagonal contribution for the auto form factor differing from (l5j) by the absence 
of the sign factor (— 1)"^ is known to stem from pairs of the so called orbit-partners. 
A partner 7' of an orbit 7 is the orbit which practically coincides with 7 everywhere 
but in the "encounters", i.e., a set of two or more orbit stretches almost parallel and 
abnormally close to each other. The lengths of 7 and 7' can be so close that their action 
difference can be comparable to or smaller than h. Consequently the corresponding 
contributions to the auto form factor will avoid the destructive interference which, in 
the case of randomly composed pairs, would lead to annihilation of the contribution. 
The leading off-diagonal correction is due to the so called Sieber-Richter pairs [T7] : 
higher order terms giving the auto form factor as a series in powers of r are calculated 
in [18] . Analytically continuing the result one gets the auto form factor for times smaller 
than the Heisenberg time, i.e., for r < 1; for larger times a different approach is needed 
(see the following subsection). 

Let us now address the cross form factor. Again, contributions of randomly 
composed pairs are expected to average to zero. Limiting summation only to pairs 
of the orbit-partners we may reduce ([5]) to 



The inner sum is taken over all partners 7' of the orbit 7 including the diagonal ones. 
Since the partner orbit coincides with the original orbit everywhere except the relatively 
short encounter stretches, the Maslov indices of 7 and 7', coincide, apart from the 
contribution of reflections from the Dirichlet boundary. Again, we assume that there 
is no correlation between parity of the number of strikes of the orbit 7 against the cut 
and the value of the sum over partners of 7 multiplied by its stability factor. Hence 
the expected value of Kg^ (r), with the off-diagonal contributions taken into account, is 
zero, at least for r < 1. 

3.4. Times larger than the Heisenberg time 

The fact that maxima (and in the case of GSE, even singularities) of form factors can 
occur at the Heisenberg time, makes it necessary to investigate the long-time behaviour 
in the discrete symmetry case. Semiclassical investigation of spectral auto and cross 
correlation beyond the Heisenberg time was started by Bogomolny and Keating|19j. 
Complete expansion of the autocorrelation functions using the formalism of the 4- 
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determinant al generating functions is given in [201 EI]; the results of the latter can 
be generalized to cross-correlations. 

The appropriate generating function of the two-point cross-correlation functions for 
the spectra of Hamiltonians H'^^\ H^'^'> is a ratio of four spectral determinants, 

^ _ / det + 5cA/27r - H^')) det [E + A/27r - H^^)) \ 
Z [SA, SB, ec, So) - ^ ^ ^^^/2vr - iJd)) det {E + sbN'^^ - 

where {■ ■ ■) ^ usual indicates averaging over the center energy E and smoothing over 
small windows of the energy offsets from E\ unlike the generating function in the problem 
of autocorrelation it is not symmetric with respect to the exchange of Sc and Sd- Using 
the identity 

oE 

and proportionality between ImTr [E — H) ^ and the level density p {E) we obtain the 
cross-correlation function as the second derivative, 

R^^'^He) =Re hm (Tr (E + ^-ifW^ \i(e-^^^ 




7?^+o \ \ 27r / V 27r 

= 2Re lim - — - — Z^^^^ (e:^, ^b, ^c, ^d) 
= e + irj . 

Here ?7 is a small positive number, and the symbol "||" denotes the substitution of the 
arguments, 

Ea = Sb = 



The semiclassical approximation for the generating function follows from that for 
the spectral determinants. Using the so called " Riemann-Siegel lookalike" asymptotics 
of det {E — H) introduced by Berry and Keating |22l |23l El] one can show that the 
semiclassical estimate of Z^^'^^ consists of two terms. One of these substituted into ffTTl) 
leads to a correlation function whose Fourier transformation reproduces the form factor 
for small times. The second term in Z^^'^^ is responsible for the behavior of the form 
factor for times larger than T^- Each term in the semiclassical generating function 
is obtained as an infinite product over the periodic orbits, averaged over the central 
energy E. The leading order contribution is obtained by the diagonal approximation 
which, in the formalism of the generating functions, means that the factors associated 
with different periodic orbits are assumed independent (with the obvious exception of 
the mutually time-reversed orbits); consequently averages of products over the periodic 
orbits are replaced by the products of averages. 
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The diagonal approximation for the generating function describing the 
correlation is deduced in the same way, step by step, as in the case of autocor- 
relation; for the lack of space we address the readers to the papers [201 El] • At a certain 
stage one gets. 



Zi = ei(^^-=s-£c+£i5) JJ exp g^j „^ 

7 

Zii = e^(^^-=s-£D+ec) JJ exp g^jj .^ . 

7 

The factor g in the exponents is equal to 2 for systems with time reversal allowed (like 
the cardioid billiard) where it reflects the existence of pairs of mutually time reversed 
orbits; in the absence of time reversal g = 1- The exponents are given by 

^jj,^ = - \F^f [e^^^^A - (-1)"^ e*"^""] [e"^"^"^ - (-1)''^ e'"-<'^] . 

We denote here by = T^/Tu the period of a periodic orbit 7 in terms of the Heisenberg 
time; is its stability coefficient. As before is the number of visits of the orbit 7 
to the cut in the desymmetrized billiard; if the factors (—1)"^ were dropped we would 
return to the case of autocorrelation. 

It remains to perform summation over the periodic orbits in the exponents of 
Zj, Zii. Again, we assume absence of correlation between parity of n^, on the one 
hand, and the orbit period and the stability coefficient, on the other hand. Then the 
summands containing the factor (—1)'^^ disappear after summation over 7 such that 



7 

^^ua ^ - XI l^-^l' [e'"-^"^"'''^ + e'"-("°-^^)] , (11) 
7 7 

consequently the term Zj in the generating function is an irrelevant constant. 

The sum in the second line of (ITTl) can be evaluated using the sum rule of Hannay 
and Ozorio di Almeida [25] for the stability coefficients F^, 

2_ AT 

7I 



r<r^<T+Ar 

which gives 



T 



^ In [{ea - Be) {bd - es)] + 2 In tq + 27 - ivr 
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with 7 denoting Euler constant. The lower integration hmit is tq = Tq/Th with Tq 
standing for some minimal period above which periodic orbits behave approximately 
ergodically; in the semiclassical limit tq tends to zero. The second term in the generating 
function is thus, 

Zn SB, ec, Sd) oc et(^^-^^-^^+^^) [r^ {sj, - sc) [ed - e^)] ' . 

The respective two-point correlation function is zero in the presence of time reversal 
when g = 2. Without time reversal {g = 1) the term Zjj contributes like const TqC*^^ to 
the cross-correlation function or, equivalently, like a delta-peak to the time dependent 
cross form factor Kg^ (r) at r = 1, i.e., at the Heisenberg time. However, since the term 
is proportional to Tq this contribution vanishes in the semiclassical limit. Therefore, as 
could be expected, correlations between the even and odd part of the spectrum do not 
exist for times either smaller or larger than the Heisenberg time. 



3.5. General discrete symmetries 

In this section we give a formal treatment of the problem in the presence of more 
complicated discrete symmetry groups G. The level density p (E) falls into a sum of 
subdensities pm (E) connected with different irreducible representations (IR) m of G. 
The semiclassical approximation for the fluctuating part of pm {E) is provided by the 
symmetry-adapted Gutzwiller trace formula ^U\. In the following, we shall follow the 
ideas expounded in [TO], so that our treatment includes the case in which the reduction 
of the configuration space M of the system to a fundamental domain M, which tesselates 
the full M when the symmetry operations are applied to it, is difficult or impossible. 
Let us look at the cross form factor defined as in ([2]) 

k,l 

throughout this section the argument of the form factor r = t/Tn, i. e., time measured 
in units of the Heisenberg time. What we want to show is that this quantity vanishes. 
At the abstract level this can be rewritten as 

= [Tr (P^,e-^^*/^'P^J]^ ■ [Tr (P„,e-^^*/''P„,)] (13) 

where P^ is the projector of the Hilbert space on the IR m. It is given by the expression 

Pm = T^\Y.^"^^3)9 (14) 

' ' geG 

where Xm{g) is the character corresponding to the IR m and g denotes the group element, 
which is assumed to act directly on the elements of the Hilbert space. 

As is well known, traces can be estimated semiclassically as integrals over coherent 
states \p, q). Assuming that the symmetry operations g have a well defined semiclassical 
limit, we know how they act on the phase space as well as on the Hilbert space, and 
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thus are able to define such quantities as g{p, q). One therefore obtains 

X XmA9'){P2,q2\e-''''/'\g'{p2,q2)) (15) 
where we have used the identity Tt{PAP) = Tt{AP) where P is an arbitrary projector. 
Since the \p, q) are well-localized in phase space, the kind of expressions occurring in 
f lTSj) only differs from zero when the phase space point (p, q) is (at least approximately) 
carried over to g{p, q) in time t by the classical evolution. In any case, this reasoning is 
legitimate when t is less than the Ehrenfest time. Since, in this section, we shall limit 
ourselves to the diagonal approximation, which is only valid when times are not too 
large, this is a reasonable approach. Now we must account correctly for all the various 
contributions to this integral. 

To this end, let us start by looking at the quantity Tr (Pmie~*^*''''Pmi) • Let us 
define as a generalized periodic orbit of period t any segment which connects a phase 
space point (p, g) with one of its symmetric partners g{p,q) over time t. We call such 
orbits 5f-orbits, making explicit reference to the corresponding group element. Let A be 
a (yf-orbit. One then additionally defines gx to be the minimal element of the symmetry 
group G such that gx{p,q) lies on the orbit A going through {p,q) and such that no 
other element g' of G has the property that g'{p, q) lies on the orbit A between (p, q) 
and g\{p,q). Note that in general g\ ^ g. There is, however, as shown in [10], always 
an integer n such that g = g^. 

Let (p, q) lie on a (^-orbit A. One then has 

l,p,q\e-^'''l^\g{p,q)) = A^e^'-l^ (16) 

where Ax is related to the stability of the orbit in the usual way and the phase e'"^^/^ 
contains all the information on the action of the (7-orbit as well as its Maslov phases. 
Integrating the l.h.s. of f|T6|) over p, q leads to 

1 t °° 

Tr (P„,e-^*/^^P^,) = ^ E 7^ E E ^>^^'''"'x.nM)h.i (17) 

I I geG ^ ' n=l A 

where n is the number of iterations of the primitive element leading to g. The factor 
t arises from the fact that an integration over the whole segment from (p, g) to g{p,q) 
is possible and yields the same result as starting with {p,q). The sum still runs 
independently over the various symmetric "copies" of the (/-orbit A, in particular over 
its iterations which lead eventually to a bona fide periodic orbit. Neglecting repetitions, 
that is n 7^ 1, yields 

Tr (P„,e-^^*/^P™,) = ^^^^ E ^Ae'^^/^Xn.. (^?a) • (18) 

Let us now multiply this with the complex conjugate of the corresponding expression 
for mi. This gives: 

K^-^-'^^t) = (^^^^^E^^^>^^'''"''^^^'^-i (9,)Xm, (gx) . (19) 
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One now introduces the diagonal approximation, which involves keeping in this sum 
only such terms as have a phase equal to one, that is, only to sum over these pairs of A 
and fi such that S\ = S^j. Since we consider the different iterations of a g orbit to be 
different, and since we have already discarded orbits which are non-primitive, and since, 
in the same vein, we discard such orbits as have a non-trivial symmetry (that is, we 
assume that the only h E G such that h{p, q) = (p, q) is the identity), it is clear that the 
total number of such pairs is [Gp independently of A and fi: there are |G| degenerate 
A's and the same number of degenerate /I's. The final result is hence 

^^'"^•™^Hr) = ^ Yl l^^l'^-i (9x) X™. (gx) (20) 

A 

where summation is carried out now over one representative (7-orbit A among the \G\ 
degenerate ones. In other terms, the sum corresponds to a desymmetrized problem. 
Keeping all ^f-orbits would require a compensating factor of in front of the r.h.s. 

of 

One now rewrites this as 

^ ^ a A 
and we now make a final assumption altogether similar to those that we have been 
making throughout in this paper: namely, we assume that when the (7-orbits are 
sufficiently long, they are equidistributed on all elements of the group G. This leads to 
stating that the sum over A is in fact independent of g. This is, of course, sufficient to 
guarantee that K^"^'^'"^'^\t) = whenever mi 7^ m2. 



4. Continuous symmetry 

In the presence of continuous symmetry, the problem has one or several integrals of 
motion due to Noether's theorem. We focus on one of them which can be chosen as 
an action and whose eigenvalues in the semiclassical limit are given by hm, where m is 
integer. Our question is whether correlation may exist between the energy spectra of 
two operators, Hm and Hm+i- Here we denote by Hm the Hamiltonian H reduced to 
the appropriate symmetry subspace. In the following, we limit ourselves to the case in 
which the dynamics of these Hamiltonians are all overwhelmingly chaotic in the energy 
range of interest. 

A simple example which we shall use to illustrate our argument is a system with 
axial symmetry in a 3d space. Introducing cylindric coordinates we can separate 
the angular part of the eigenf unctions writing \l/ = e*™'^?/' (p, z) where if) (p, z) is an 
eigenfunction of the Hamiltonian Hm (p, z) in the subspace of states with the given z- 
component of the angular momentum = hm. As stated above, we assume that the 
2d classical motion described by this Hamiltonian is chaotic. A moment's reflection will 
show that this example contains all the relevant features of the more general situation 
described at the beginning of this section. 
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In this case, the Hamiltonians Hm and Hm+i correspond to the radial Hamiltonians 
with magnetic quantum number m and m + 1 respectively. The magnetic quantum 
number enters these Hamiltonians as a parameter, and the well-known theory of 
parametric correlations can be applied. Special treatment is necessary when m = 
or m ~ 1. The difference between the two Hamiltonians due to the centrifugal term 
W = Ll/2Mp^ (M is the mass of the particle) is then of the formal order f^. However the 
energy shifts associated with the change of m by 1 are of the order h\ this discrepancy 
is due to singularity of the centrifugal term at p = 0. The appropriate Gutzwiller 
expansions in this case use the set of classical periodic orbits calculated for = 0, 
i.e., m- independent. The m-dependence of the spectra is taken into account via an 
appropriate phase added to the Maslov phase; the phase increment occurs each time the 
orbit hits the axis p = 0, see the reviews |26l [27]. Since this additional phase is pseudo- 
random, the cross-correlation functions between the spectra with different m ~ 1 vanish 
due to the same mechanism as between the spectra of different parity in the discrete 
case. 

Let us now turn to the case of \m\ ^ 1; more precisely consider the semiclassical 
limit /z — )■ 0, m — i- oo such that the angular momentum projection Lz = hm remains 
constant. In this case the centrifugal barrier W would not let the particle approach 
the z -axis; the singularity at p = will thus be deep in the classical shadow. The 
change of m by 1 brings about change of the Hamiltonian of the order h, namely 
by AHm = AW = hdW/dLz- Before comparing the eigenvalues at different m one 
must subtract the overall drift of the spectrum upwards caused by the increase of the 
centrifugal energy W. Note that this issue is similar to the one caused by the short 
periodic orbits in the discrete case: it arises from a difference in the one-particle densities 
of the two systems, and is an issue which can also be settled by means of an appropriate 
unfolding procedure. This quantity can be evaluated semiclassically as the derivative of 
its average over the energy shell. 



where dT is the element of the phase space. Denoting by 5*0 a typical classical action 
within the given (highly excited) energy range, we can introduce a dimensionless small 
parameter ^ = H/Sq and write 



consider the general problem of parametric correlation in the family of Hamiltonians 
= Hq+^H' where ^ is small and dimensionless; the Hamiltonians Ho, H' are assumed 
to have finite classical counterparts. Let us consider the cross form factor between the 
spectra of Hm and H^+i in the semiclassical domain. We shall assume that ^ is so small 



{W) 



J dT6 {E - Hm) W 
f dT6 {E - Hm) 




To appreciate the impact of the perturbation associated with Hm H, 



m+l 
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that the periodic orbits of Hm, H^+i practically coincide; only their action difference 
has to be taken into account since it is referred to h. 

The additional classical action due to the perturbation H' over a stretch of the 
orbit with duration from ti to ^2 would be 

AS~e r H'dt. 

Now the well-established approach is to model the accumulation of AS on a trajectory of 
a chaotic system by a Gaussian stochastic process [28]; below ((...)) denotes averaging 
over such a process . This can be understood if we mentally divide the periodic 
orbit into pieces of duration larger than the Lyapunov time = where A is the 
Lyapunov constant; the action differences accumulated at two such stretches would be 
uncorrelated. Let t = t2 — ti ^ and assume that there is no systematic action 
increment, {{AS)) = 0. The averaged second moment of AS* will grow proportional to 
t , 

{{AS')) = aet 

where a is a purely classical, system- specific, parameter; its dimensionality is 
squared action over time. Remembering that in a Gaussian process ((exp(i$))) = 
exp (— (($^)) /2) we obtain the average of the exponential of the phase difference 
accumulated on a periodic orbit with the period T^, 

Now let us use the semiclassical expression (jlj) for the cross-correlation form factor 
which we shall now denote (r) . Limiting ourselves to the diagonal approximation we 
would get, for the time T = tTh, 

K,,..^Ar)-\ E (22) 

T<T^,<T+S 

= e-'^fi^i^diag (r) . (23) 

The quantity -f^diag (t) is the autocorrelation form factor in the diagonal approximation, 
i.e., 2r or r depending on the presence or absence of the time reversal. The perturbation 
^H' leads thus to the decay of the form factor with the characteristic time Tjecay = 

Applying this result to the correlation between the spectra belonging to two 
neighboring values of the quantum number m when = h/So we obtain the decay 
time 

decay 

a 

This time is classical, i.e., finite in the semiclassical limit, however it could be much 
larger than the ballistic time scale, i.e., the period of the shortest periodic orbits (time 
on which correlations are expected anyway). 
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The corresponding dimensionless Tdecay = ^decay/^H turns into zero in in the hmit 
— > 0. We therefore find that the cross-correlation between spectra corresponding to 
different symmetries vanishes in the universal regime. Considering that the periodic 
orbit expansions of the type (jl]) — and above all the averaging procedure over large 
numbers of periodic orbits having nearly the same action — are justified only for times 
exceeding the Ehrenfest time Te ~ In Sq/H, the vanishing of the cross-correlation 
only holds for times larger than Te- For times shorter than Te, and hence in particular 
for times of the order of classical times, correlations are expected to arise, based on the 
fact that both symmetry sectors have the same periodic orbits. Such correlations will 
depend on the specific features of the periodic orbits and will thus be highly system- 
dependent. A possible connection with a recent discussion of the transition between 
short time effects and the universal regime [29] is left for future investigation. 

5. Conclusions 

In this paper we have given a justification to the common assumption that spectra 
corresponding to different irreducible representations of a symmetry group of a given 
quantum chaotic problem are independent in the semi-classical limit. More precisely 
we have shown in the semi-classical limit that the two-point cross correlation vanishes, 
a property also sometimes called weak independence. For the specific case of a mirror 
symmetry we showed that the cross correlation function is zero not only in the diagonal 
approximation and for short times. The non diagonal terms of the trace formula were 
also evaluated and the semi-classical form of the spectral determinant was used to show 
that our result holds also beyond the Heisenberg time. For arbitrary discrete symmetries 
and for continuous symmetries we limited our argument to the diagonal approximation, 
but the arguments given for mirror symmetries carry over. 

An potentially interesting line for future work appears in continuous symmetries, 
where the trace formulae reveal an unusual correlation on a classical time scale. The 
approximation used is clearly not valid on this time scale, but the result might still be a 
hint, that such correlations exist and are different from the obvious ones resulting from 
short orbits. 
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